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We study a contact process running in a random environment in 
Z'' where sites flip, independently of each other, between blocking and 
nonblocking states, and the contact process is restricted to live in the 
space given by nonblocked sites. We give a partial description of the 
phase diagram of the process, showing in particular that, depending 
on the flip rates of the environment, survival of the contact process 
may or may not be possible for large values of the birth rate. We 
prove block conditions for the process that parallel the ones for the 
ordinary contact process and use these to conclude that the critical 
process dies out and that the complete convergence theorem holds in 
the supercritical case. 

1. Introduction. We consider the following version of a contact process 
running in a dynamic random environment in Z*^. The state of the process is 
represented by some i] & X = { — 1,0,1} , where sites in state are regarded 
as vacant, sites in state 1 as occupied and sites in state —1 as blocked (that 
is, no births of I's are allowed on that site). The process Tjt is defined by the 
following transition rates: 

— > 1 at rate (3fi 

1 — > at rate 1 
0, 1 — > —1 at rate a 
— 1 — > at rate a6 

where /i is the fraction of occupied neighbors at distance 1. 

In words, the —I's define a random environment in which each site be- 
comes blocked at rate a and flips back to being unblocked at rate a5, while 
the I's behave like a nearest neighbor contact process with birth rate (3 in 
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the space left unblocked by the environment. Observe that when an occupied 
site becomes blocked, the particle is killed. This version is simpler than the 
alternative in which only O's can turn to — I's (mainly because our process 
satisfies a self-duality relation, see Proposition 2.2). However, we feel that 
our choice is natural: If a site becomes uninhabitable, the particles living 
there will soon die. 

Ever since it was introduced in Harris (1974), the contact process has 
been object of intensive study, and many extensions and modifications of the 
process have been considered. In particular, the literature includes several 
different versions of contact processes in random environments. One class of 
these processes corresponds to contact processes where the birth and death 
rates are not homogeneous in space, and they are chosen according to some 
probability distribution, independently across sites, and remain fixed in time 
[see, e.g., Bramson, Durrett and Schonmann (1991), Liggett (1992), Andjel 
(1992) and Klein (1994)]. The main question for this class of processes is to 
determine conditions on the parameters that guarantee or preclude survival. 

A different class of models, which are somehow closer to the process we 
consider, have two species with different parameters or ranges, but one of 
them behaves independently of the other while the second is restricted to live 
in the space left by the first. These processes were studied in Durrett and Swindle 
(1991), Durrett and M0ller (1991) and Durrett and Schinazi (1993). The re- 
sults in these papers (mainly bounds on critical parameters for coexistence 
and complete convergence theorems) are asymptotic, in the sense that they 
are proved when the range of one or both types is sufficiently large. 

The process we consider differs from both of the types of examples men- 
tioned above: The random environment is dynamic and it behaves indepen- 
dently across sites. An example of a spin system running in this type of 
environment was studied in Luo (1992), and corresponds to the Richardson 
model which would result from ignoring transitions from 1 to in our pro- 
cess. Another example was studied recently in Broman (2007), where the 
author considers a process in which the environment changes the death rate 
of the contact process instead of blocking sites. The dynamics of the process 
^Ifsf introduced there are the same as those of our process if (5i = oo. The 
author considers this tool in the study of his process, but the re- 

sults of the paper focus on the case 5i < oo. We will use one of his results to 
give a bound on a part of the phase diagram of our process in Theorem 1. 

As mentioned above, the — I's evolve independently of the I's. They follow 
an "independent fiip process" whose equilibrium is given by the product 
measure 





This process is reversible, and its reversible measure is given by fip. 
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In Section 2.1, we will construct our process using the so-called graphical 
representation. A direct consequence of the construction will be that rjt 
satisfies some monotonicity properties analogous to those of the contact 
process. (Here and in the rest of the paper, when we refer to the contact 
process we mean the "ordinary" nearest-neighbor contact process in Z'^.) 
We consider the following partial order on configurations: 

(1.1) V^<v'^ ^ ri^{x)<rf'{x) Vx G Z'^. 

With this order, our process has the following property: Given two initial 
states tjq <r]Q, it is possible to couple two copies of the process rj^ and ry| 
with these initial conditions in such a way that r]} < rjf for all t > 0. We 
will refer to this property as attractiveness by analogy with the case of spin 
systems [this property is sometimes termed monotonicity, see Sections II. 2 
and III. 2 of Liggett (1985) for a discussion of general monotone processes 
and of attractive spin systems, resp.]. 

For A C Z*^, we define the following probability measure ua on X: — I's 
are chosen first according to their equilibrium measure fj,p and then I's are 
placed at every site in A that is not blocked by a —1. These measures are 
the initial conditions for r]t that are suitable for duality. 

Let z^ = z^0, which corresponds to having the —I's at equilibrium and 
no I's. Let also V be the limit distribution of the process when starting at 
the configuration having all sites at state 1, which is obviously the largest 
configuration in the partial order (1.1). We will show in Proposition 2.1 
that this limit is well defined and it is stationary, and that and V are, 
respectively, the lower and upper invariant measure of the process (i.e., the 
smallest and largest stationary distribution of the process) . 

We will say that the process survives if there is an invariant measure u 
such that 

z^({r7 : ri{x) = 1 for some x G Z'^}) > 0, 

or, equivalently, if z^t^I^ (we remark that, as a consequence of Theorem 2, 
every invariant measure for the process is translation invariant, so the above 
probability is actually 1 whenever it is positive). Otherwise, we will say that 
the process dies out. We will see in Section 4 that this definition of survival 
is equivalent to the following condition: The process started with a single 1 
at the origin and everything else at —1 contains I's at all times with positive 
probability. 

A second monotonicity property that will follow from the construction of 
r]t is monotonicity with respect to the parameters P and 6: 

(i) If a and 5 are fixed, and for some /? > the process survives, then 
the process also survives for any /?' > (3. 
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(ii) If a and /? are fixed, and for some 5 > the process survives, then 
the process also survives for any 5' > 5. 

These properties follow easily from standard coupling arguments. We will 
denote by Pc = Pc{oi, 6) € [0, oo] the parameter value such that, fixing these a 
and 6, r]t survives for P > (3c and dies out for (3 < f3c- We define dc = dc{a, 13) 
analogously. 

Our first result provides some bounds on the critical parameters for sur- 
vival. Let be the critical value of the contact process in [here we are 
taking the birth rate (3 to be the total birth rate from each site, so each site 
sends births to each given neighbor at rate (3/{2d)\. 

Theorem 1. 

(a) If (3 < {a + l)/3c^, then the process dies out. 

(b) There exists a 5p > such that for any 5 < 6p the process dies out 
(regardless of a and j3). 

(c) Let 



If A(q!,/3,(^) > (a + 1)13^ , then the process survives. 

Part (a) of the theorem is trivial because the I's die at rate a + 1. For 
part (b), observe that if the complement of the set of sites at state — 1 
does not space-time percolate, then each 1 in the process is doomed to 
live in a finite space-time region, and then the process cannot have I's at 
all times when started with finitely many occupied sites. We will show by 
adapting arguments in Meester and Roy (1996) that, with probability 1, no 
such space-time percolation occurs if 5 is small enough. For part (c), we will 
use Broman's result to obtain a suitable coupling with a contact process 
with birth rate A(a,/3,(5) and death rate a + 1. 

In particular. Theorem 1 implies that if 5 is large enough then /?c(a,(5) < 
oo, and in fact 5 > ^^Pc^ is enough. To see this, observe that 



whenever the above condition on 6 holds. Then part (c) of the theorem 
implies that the process survives for these choices of a and 5 and large 
enough (3. Another consequence is that 6p < Indeed, ii 6 > then 
6 > for large enough a, and the previous property implies that the 

process survives for these choices of a and 6 and large enough p. 

A significant difficulty in giving a more complete picture of the phase 
diagram of r]t is that we lack a result about monotonicity with respect to 




lim A(a, P, 6) = a5 > {a + 



[Cp 
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a analogous to the properties (i) and (ii) (monotonicity with respect to /? 
and S) mentioned above. Observe that the equihbrium density of nonblocked 
sites is independent of a, but the environment changes more quickly as a 
increases. Simulations suggest that if /? and b are given and the process 
dies out at some parameter value a, then it also dies out for any parameter 
value ol > a (note that part (a) of Theorem 1 says that the process dies 
out at least for all a large enough). But the usual simple arguments based 
on coupling do not work in this case, since increasing a increases both the 
rate at which sites are blocked, which plays against survival, and the rate 
at which sites are unblocked, which plays in favor of survival, and we have 
not been able to find an alternative proof. 

The second part of our study of r\t investigates the convergence of the 
process and the structure of its limit distributions. For t] X, we will write 
rj = {A, B), where 

A = {xeZ'^:r]{x) = l} and B = {x eZ'^ :r]{x) = -1}. 

7]^ = (A^,B^) will denote the process with initial distribution fi, and we 
will refer to B^" (or Bt if no initial distribution is prescribed) as the envi- 
ronment process. Observe that the dynamics of the environment process are 
independent of the I's in r]t. 

Theorem 2. Denote by t = inf{t > : At = 0} the extinction time of 
the process. Then for every initial distribution fi, 

rji^ =^ P^(r < oo)u + P^(r = oop, 

where the limit is in the topology of weak convergence of probability measures. 

This result, which is usually called a complete convergence theorem, im- 
plies that all limit distributions are convex combinations of and V. Thus, 
the only interesting nontrivial stationary distribution is V. 

The proof of Theorem 2 relies on extending for rjt the classical block con- 
struction for the contact process introduced in Bezuidenhout and Grimmett 
(1990), so that we are able to use the proof of complete convergence for the 
contact process to prove the corresponding convergence of the contact pro- 
cess part of rjt- As a consequence of this construction, we will obtain, just 
as for the contact process, the fact that the process dies out at the critical 
parameters /3c and 6c (see Corollary 4.4). The arguments involved in this 
part will depend heavily on a duality relation which will be developed in 
Section 2.2. 

The rest of the paper is devoted to the proofs of the two theorems. Sec- 
tion 2 describes the construction of r]t and presents some basic preliminary 
results. Theorem 1 is proved in Section 3. In Section 4 we obtain the block 
conditions for the survival of the process. Finally, in Section 5 we use duality 
and the conditions obtained in Section 4 to prove Theorem 2. 
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2. Preliminaries. 

2.1. Graphical representation and monotonicity. The graphical represen- 
tation is one of the basic and most useful tools in the study of the contact 
process and other interacting particle systems. It will allow us to construct 
our process from a collection of independent Poisson processes and obtain a 
single probability space in which copies of the process with arbitrary initial 
states can be coupled. We will give a rather informal description of this con- 
struction, which can be made precise by adapting the arguments of Harris 
(1972). We refer the reader to Section III. 6 of Liggett (1985) for more details 
on this construction in the case of an additive spin system. 

The construction is done by placing symbols in Z'^ x [0, oo) to represent 
the different events in the process. For each ordered pair x, y E Z'^ at distance 
1, let N^'^ be a Poisson process with rate P/{2d), and take the processes 
assigned to different pairs to be independent. At each event time t of N^'^ , 
draw an arrow in Z'^ x [0,oo) from {x,t) to {y,t) to indicate the birth 
of a 1 sent from x to y (which will only take place if x is occupied and 
y is vacant at time t). Similarly, define a family of independent Poisson 
processes {U^'^)xeZ'i with rate 1 and for each event time t of U^'^ place a 
symbol *i at {x,t) to indicate that a 1 flips to (i.e., that a particle dies). 
To represent the environment, consider two families of independent Poisson 
processes {V^)rc^^d and (C^~^'^)a;ez<^ with rates a and a6, respectively. For 
each event time t of V^, place a symbol •_! at {x,t) to indicate the birth of 
a —1 (i.e., the blocking of a site) and for each event time t of place 
a symbol *_i to indicate that a —1 flips to (i.e., the unblocking of a site). 

We construct rjt from this percolation structure in the following way. Con- 
sider a deterministic initial condition rjQ and define the environment process 
Bf by setting 'i]t{x) = —1 when {x,t) lies between symbols •_! and *_i (in 
that order) in the time line {x} x [0, oo), and also if r]Q{x) = — 1 and there 
is no symbol *_i in that time line before time t. Having defined Bt, we say 
that there is an active path between (x,s) and {y,t) if there is a connected 
oriented path, moving along the time lines in the increasing direction of time 

and passing along arrows which crosses neither symbols *i nor space- 
time points that were set to —1. The collection of active paths corresponds 
to the possible space-time paths along which I's can move, so we define At 
by 

At = {y 1^'^ :3x (z Aq with an active path from (x, 0) to {y, t)}. 

The arguments of Harris (1972) imply that this construction gives a well- 
defined Markov process with the right transition rates. Moreover, the same 
realization of this graphical representation can be used for different initial 
conditions, and this gives the coupling mentioned above [see Section III. 6 
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in Liggett (1985) for more details on this coupling in the case of a spin 
system]. For the rest of the paper, we will implicitly use this "canonical" 
coupling every time we couple copies of rjt with different initial conditions. 
The attractiveness property mentioned in the Introduction follows directly 
from this construction, and the monotonicity properties with respect to /? 
and 5 can be obtained by a simple modification of this coupling (analogous 
to what is done for the contact process). 

Recall the definition of the partial order on configurations given in (1.1). 
Clearly, 

'n^<'n^ ^ Q A"^ and ^ B"^ . 

For probability measures on X, which we endow with the product topology, 
we consider the usual ordering: fii < ^2 if and only if / / d^ii < J f d^2 for 
every continuous increasing / : X — > M. We recall that the property ^1 < ^2 
is equivalent to the existence of a probability space in which a pair of random 
variables Xi and X2 with distributions ^1 and 112 can be coupled in such a 
way that Xi < X2 almost surely [see Theorem II. 2. 4 in Liggett (1985)]. We 
will use this fact repeatedly, and for simplicity we will say that X2 dominates 
Xi when this condition holds. We will also use this term to compare two 
processes, so saying that r]^ dominates r]^ will mean that the two processes 
can be constructed in a single probability space in such a way that 77^ < r]^ 
for ah t > 0. 

The attractiveness property allows us to obtain the lower and upper in- 
variant measure of the process. 

Proposition 2.1. Let xz^ be the probability distribution on X assigning 
mass 1 to the all 1 's configuration, and let S{t) be the semigroup associated 
to the process. Define 

F= lim XziS{t), 

t— >oo 

where the limit is in the topology of weak convergence of probability mea- 
sures. Then V is the upper invariant measure of the process, that is, V is 
invariant and every other invariant measure is stochastically smaller than 
V. Moreover, 

V = lim UvdSU). 

t^oo 

Analogously, 
is the lower invariant measure of the process. 
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Proof. Since /ip is invariant for the environment and the empty state 
is a trap for the I's, v_ is invariant. It is the lower invariant measure because 
every invariant measure has fXp as its projection onto the environment, and 
Vgj is the smahest probabihty measure on X having /ip as its marginal on 
the -I's. 

For the other part, standard arguments imply that the limit defining 
u exists and is invariant [see, e.g., Sections I.l and III. 2 in Liggett (1985)]. 
Since Xz'* is larger than any other measure on X follows by attractiveness 
that V is the largest invariant measure. 

Now let V* = \\mt^aoVj^dS{t). As above, v* is well defined and invariant, 
so to prove that v* =V it is enough to prove that v* is larger than any other 
invariant measure. If v is any invariant measure, its projection onto the —I's 
must be /ip, so for any continuous increasing /, 

j /dz. = E'^(/(7?o))=E-(/(r?t)) 

<W^'^{f{r^t))= j fd[u^,S{t)]^-^ j fdu*. □ 

2.2. Duality. The dual process (r}*)o<s<f = (^s, -Bs)o<s<t is constructed 
using the same graphical representation we used for constructing r]t. Our 
duality relation will require that the process be started with the environment 
at equilibrium. The dual processes will also be started with measures of the 
form uc, for C C Z*^, and the dual process started with this distribution will 
be denoted by (C^'*)o<s<t- 

Fix t > 0, and start by choosing Bq according to /_fp. Then run the envi- 
ronment process forward in time until t, using the graphical representation. 
This defines (-Bs)o<s<t- The dual environment is given by = Bt^s- Now 
place a 1 at time t at every site x G C \ Bq, that is, every site in C which is 
not blocked by the environment at time t. This defines Aq'^'^, and by the sta- 
tionarity of the environment process we get an initial condition {Aq'^'^,Bq) 
for the dual chosen according to i^c- Having defined Aq*^'* and (-B*)o<s<t, 
we define the 1-dual by 

i^c,t = {y^Z'^:3xe io'^'* with an active path from {y, t - s) to (x, t)}. 

That is, the 1-dual is defined by running the contact process for the I's 
backwards in time and with the direction of the arrows reversed. An active 
path in rjt from (y, t — s) to (x, t) will be called a dual active path from (x, t) 
to {y, t — s) in the dual process. 

We could have defined the dual by simply choosing a random configuration 
at time t according to vc and then running the whole process backward. The 
idea of the preceding construction is to allow coupling the process and its 
dual in the same graphical representation in such a way that the initial state 
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of the environment for rjs is the same as the final state of the environment for 
57* (that is, Bq = Bj). This allows us to obtain the following duality result: 

Proposition 2.2. For any A,C,D C 1/^, 

(2.1) P''^ {AtnC ^0,Btr\D^0) = P^c^ {A\r\A^0,Blr\D^0). 

Moreover, r]t satisfies the following self-duality relation: if A or C is finite, 
then 

(2.2) F''-^{AtnC^0,Btr\D^0)=F'''^{Atr\A^0,Bor\D^0). 

Proof. The first equality follows directly from coupling the process and 
its dual using the same realization of the graphical representation. Indeed, 
if we use this coupling then, by definition, 

r{B^^c;t = B'^A^ fQj. g^g^y < s < i) = 1. 

Calling £ the cj-algebra generated by the environment process, observe that 
our construction implies that 

W''^{Atr\C ^ 0\£) = ¥''c {A\f^ A^ 0\£). 

Therefore, 

P'"^(At n c / 0, n / 0) = ¥r^{¥{At n c 0\£),Bt n / 0) 

= {¥{A\ r\A^0\£),B^Qr\D^0) 

= ¥'"^{Alr\A^0,B^QnD^0). 

Equation (2.2) is obtained from (2.1), the self-duality of the contact pro- 
cess, and the reversibility of the environment. □ 

Taking A finite and C = D = j/ in (2.2) and using the monotonicity of 
the event {At 7^ 0} in t we obtain the following: 

F^^iAt / Vi > 0) = V{{{E, F):EnA^ 0}). 

Since z7 is translation invariant, the right-hand side of this equality is positive 
if and only if A 7^ and r]t survives, that is, V ^ y_. As a consequence, 
we deduce that the following condition is equivalent to the survival of the 
process: 

For any (or, equivalently, some) finite ^4 C Z'^ with A^ 0, the 
(SI) process started at va contains I's for every t > with positive 
probability. 
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2.3. Positive correlations. A second property that is central to the study 
of the contact process is positive correlations. Recall that a probability mea- 
sure has positive correlations if for every /, g increasing, 

(2.3) J fgdf,> J fdii J gdfi. 

In the following lemma we prove a version of positive correlations for r]^^ 
with respect to cylinder functions. 

Lemma 2.3. Let f,g be increasing real-valued functions on X depending 
on finitely many coordinates. Then if fit denotes the distribution ofr^^, (2.3) 
holds with fi = nt, that is, 

(2.4) E''^firjt)g{iit))>K'Hf{r]t)WH9ir]t))- 

The same inequality holds ifvA is replaced by any deterministic initial con- 
dition. 

Proof. Since / and g depend on finitely many coordinates and every 
jump in our process is between states which are comparable in the partial 
order (1.1), a result of Harris [see Theorem II. 2. 14 in Liggett (1985)] and 
attractiveness imply that it is enough to show that the initial distribution of 
the process has positive correlations in the sense of the lemma. The result 
with UA replaced by a deterministic initial condition readily follows. 

To show that va is positively correlated, consider the process q defined 
in A" by = 1 and independent transitions at each site given by 



— > —1 at rate p 

- 1 — > at rate 1 - 

1 — > —1 at rate p 
-1 — > 1 at rate 1 - 



for X ^ A, 



for X A. 



It is clear that q converges weakly to the measure i^a- Since the initial 
distribution of q has positive correlations (because it is deterministic), (2.3) 
holds for its limit z/^, using again Harris' result. □ 

3. Survival and extinction. In this section, we prove Theorem 1. Through- 
out the proof, we will implicitly use (SI) to characterize survival. We start 
with the easy part. 

Proof of Theorem 1, part (a). Consider the process rjt defined by 
the following transition rates: 

0,-1 — > 1 at rate /3/i 

1 — > at rate 1 
0, 1 — > —1 at rate a 
— 1 — > at rate a6. 
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This process corresponds to modifying rjt by ignoring the effect of blocked 
sites on births. It is easy to couple rjt and r]t using the graphical representa- 
tion in such a way that if the initial states are the same, rjt < rjt for all t > 0. 
Therefore, it is enough to show that rjt dies out, and this follows directly 
from the hypothesis because the I's in rjt behave just like a contact process 
with birth rate (3 and death rate a + 1 . □ 

The proof of part (b) is more involved, and it is based on adapting the 
techniques of Boolean models in continuum percolation [see Meester and Roy 
(1996)]. 

Proof of Theorem 1, part (b). The idea is to show that when 5 
is small, the set of unblocked sites in the environment process Bt does not 
"space-time percolate" with probability 1. By this we mean that there is no 
infinite path in Z*^ x [0, oo) moving between nearest-neighbor sites in Z'^ and 
along time lines in the increasing direction of time that uses only nonblocked 
sites. The conclusion follows directly from this fact, since in that case, every 
1 will live in a finite space-time box, so it will not be able to contribute to 
the survival of the process. 

By a simple time change, we can consider the environment process as 
having transitions given by 

— 1 — > at rate q 
— > —1 at rate I — q, 

where q = 5/{l + 6) — > as 5 — > 0. We still consider this process as defined 
by the graphical representation, though now the symbols •_! and *_i appear 
at rate 1 — q and q, respectively. 

Take the percolation structure given by the graphical representation and 
draw for every symbol *_i at a space-time point {x,t) a box of base x + 
[—2/3,2/3]'^ spanning the interval in the time coordinate from t until the 
time corresponding to the next symbol (i.e., these boxes span intervals 
where the sites are not blocked). Then since the environment process is 
translation invariant, the O's will almost surely not space-time percolate if 
and only if 

(3.1) p(|>V| = oo) = 0, 

where W denotes the connected component of the union of the boxes that 
contains the origin at time 0, and |W| denotes the number of boxes that 
form this cluster. 

To prove (3.1), we compare this continuum percolation structure with a 
multitype branching process X = {Xn,i)n,ieN- The first step in the compari- 
son is to stretch all the boxes so that their heights are all integer- valued. It 
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is enough to show that (3.1) holds after this modification, since increasing 
the heights of the boxes increases the probabihty of space-time percolation 
of the unblocked sites. Assume that the origin is not blocked at time 0, and 
call io € N the (random) height of its associated box. For simplicity, assume 
further that all the neighbors of the origin are blocked at time 0, the ex- 
tension to the general case being straightforward. We start defining X by 
saying that the 0th generation has only one member, and it is of type io 
(i.e., Xqj = l|j=j(j}.). The box containing the origin at time is possibly in- 
tersected by boxes placed at the 2d neighbors of the origin, and these boxes 
will constitute the children of the initial member: we let Xi j be the number 
of boxes of height j that intersect the original box. We define the subsequent 
generations of X inductively: Xn+ij is the number of boxes of height j that 
intersect boxes of the nth generation and which have not been counted up 
to generation n — 1. Now let 

oo oo 
n=0 j=l 

and observe that every box in W is counted in X°°, so 

(3.2) \W\<X^ 

(recall that X is constructed from the stretched boxes). 

Our goal is to show that K{X°°) < oo. To achieve this, we will couple X 
with another multitype branching process Y = (!"«,«)", igN; which we define 
below. The details of this part can be adapted easily from the proof of 
Theorem 3.2 in Meester and Roy (1996), so we will only sketch the main 
ideas. Consider a box of height i based at [x — 2/3,x + 2/3]*^ x {t}, which we 
will denote hy B{x, t, i) . The boxes of height j that intersect this box must all 
have bases of the form [y — 2/3, y + 2/3] x {s} for some y at distance 1 of x and 
some s € (OV (t — j), t + The number of symbols *_i appearing in the piece 
{y} X (0 V (i — j) , t + i] of the graphical representation above a given neighbor 
2/ of x is a Poisson random variable with mean g[t + z — V (t — j)] <q[i+j], 
and each of these symbols corresponds to a box that intersects B{x,t,i). 
Since the probability that any one of these (stretched) boxes is of height j is 
Pj = ¥{Z S (j — 1, j]), where Z is an exponential random variable with rate 
1 — q, we deduce that the number of children of B{x,t,i) of height j is a 
Poisson random variable with mean bounded by 

(3.3) 2dpjq[i + j] < Mqijpj, 

where we used the fact that i + j < 2ij for positive integers i and j. Now 
let y be a multitype branching process where the number of children of 
type j of each individual of type i is a Poisson random variable with mean 
Adqijpj iYn^i is the number of individuals of type i in generation n). Then a 
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coupling argument and (3.3) imply that if Xo^i = Yo^j for all i > 1 then Xn^i 
is dominated by Y^^i for each n > and i > 1, and thus 

Coo oo \ 
n=Oj=l ) 

To bound this last sum, we recall a standard result in branching processes 
theory [see, e.g.. Chapter V in Athreya and Ney (1972)]: the expected num- 
ber of individuals of type j in the reth generation of Y when starting with 
one individual of type io is given by 

(3.5) ]E(y„,,-|lo,fc = l{fc=io}) = (^")io,i> 

where M is the infinite matrix indexed by N with Mjj being the expected 
number of children of type j of an individual of type i. By definition of 
Mj j = 4:dqijpj, and from this we get inductively a bound for (M")jo 

for all n > 1 , where H is a random variable with positive integer values and 
distribution given by F{H = j) = pj. Using this together with (3.4) and (3.5) 
gives 

oo / oo \ 

n=l \ j=l ) 

oo 

= 1 + 4dgioE(F) E [UqE[H^)f. 

n=0 

Observe that H is dominated by Z + 1, so E{U'^) < + 1. Hence, 

(3.7) 4«F2)<4d(^^ + ,)<l 

for sufficiently small q, and then the last sum in (3.6) converges for such q. 
This implies by (3.2) that E(|>V|) < oo, so P(|>V| = oo) = 0. □ 

Using (3.7) we can get explicit lower bounds for 6p, but these turn out to 
be rather small (around 0.02 for d = 2 and 0.01 for d = 3). 

Before proving the last part of Theorem 1, we need to introduce a result 
from Broman (2007). Let {Jt,Xt) be the process with state space {0, 1} x N 
defined as follows. Jq is a Bernoulli random variable with P(Jo = l) = l — 
P( Jo = 0) = p, and Xq = 0. The evolution of the process is given by the 
following transition rates: 

°^ *■ \ l — > at rate 7(1 -p) 
for Xf: k — > k + 1 at rate cro(l — Jt) + criJt 
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where 7,0"! > and < ctq < In words, Jt acts as the environment, start- 
ing at equihbrium and then flipping between states and 1 independently 
of Xf, while Xt is a sort of Poisson process where the rate depends on Jt. 
The next lemma recovers the part of Theorem 1.4 in Broman (2007) that 
is relevant for our purposes. We observe that the original theorem is stated 
for (To > 0, but the same proof works if do = 0. 

Lemma 3.1. Let 



Then a Poisson process Nt{a) with rate a can be coupled with {Jt,Xt) in 
such a way that if Nt{a) has an arrival at time T, then so does Xt. Moreover, 
a is the largest rate such that this coupling is possible. 

Recall that we denote 



The following result gives the coupling that we need to prove part (c) of 
Theorem 1. Its proof is very similar to that of Theorem 1.7 in Broman 
(2007); we include here a version based in the graphical representation. 

Proposition 3.2. Let S^t denote the set of occupied sites of a contact 
process with birth rate A = A(a, (3, 6) and death rate a + 1. Then the processes 
rjt and ^t can be coupled in such a way that if ^0 ^ ^O; then C At for all 



Proof. Consider the graphical representation used to construct r]t. Each 
time line defines an independent copy of the process Jt introduced above by 
identifying symbols •_! and *_i with Jt flipping to and 1, respectively, 
and setting j = a{l + 6) and p = 6/ {I + 6). Now consider the collection of 
arrows emanating from that time line ignoring arrows born at times where 
the site is blocked. By construction, this collection of arrows defines the ar- 
rival times of the process Xt associated to Jt, with ao = and ai = /3. By 
Lemma 3.1, we can construct a Poisson process Nt{X) [where A comes from 
plugging in our parameters in (3.8)] such that if this process has an arrival 
at time T, then there is an arrow at that time for rjt- 

We repeat this construction at each time line, getting an i.i.d. collection 
of Poisson processes {Nf {X)) x^i<i ■, and use this collection of processes and 
the graphical representation of rjt to construct the graphical representation 
of ^t'- for each arrival time of Nf{\) put an arrow at that time from x to the 
site pointed by the corresponding arrow in the graphical representation of 
r/t, and for each symbol *i and each symbol for rjt put a death symbol 





t>0. 
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for It is easy to see that this construction gives a graphical represen- 
tation for the desired contact process Moreover, since only the arrows 
at nonblocked sites can carry births of I's for rjt, the construction gives a 
coupling that satisfies the desired monotonicity property. These facts can be 
checked exactly as in the proof of Theorem 1.7 of Broman (2007) (there the 
processes Yf and Y"/ correspond to At and ^t). □ 

The proof of the remaining part of Theorem 1 is now straightforward. 

Proof of Theorem 1, part (c). Since ^'^"'^f^ > Ac implies that the 
contact process with birth rate X{a,f3,5) and death rate a + 1 survives, 
the coupling achieved in Proposition 3.2 gives the survival of □ 

4. Block construction. The aim of this section is to establish "block 
conditions" concerning the process in a finite space-time box that guarantee 
survival. This was first done in Bezuidenhout and Grimmett (1990). Here we 
will follow closely Section 1.2 of Liggett (1999), together with the corrections 
to the book that can be found in the author's website. 

Before getting started with the block construction, we need to obtain the 
equivalent condition for survival mentioned in the Introduction, which says 
that r]t survives if and only if the following condition holds: 

The process started with a single 1 at the origin and ev- 
(S2) erything else at —1 contains I's at all times with positive 
probability. 

The sufficiency of this condition is a consequence of (SI) and attractiveness. 
The necessity will be a consequence of the following stronger result, which is 
precisely what we will need in the proof of Lemma 4.2 below. Let XA denote 
the probability measure on X that assigns mass 1 to the configuration r] 
with tjIa = 1, vIa'^ = -1- 

Lemma 4.1. Suppose that the process survives. Then for any a >0 there 
is a positive integer n such that 

pX[-„,„]d {At^0 \/t>0)>l-a^. 

To obtain (82) from this result observe that the process started with a 
single 1 at the origin has [—n, n]"^ fully occupied by time 1 with some positive 
probability, so we can use the strong Markov property and attractiveness to 
restart the process at time 1 starting from X[-n,n]'' obtain ¥^i°'^{At ^ 
Vt > 0) > 0. Observe that the lemma is a simple consequence of duality 
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when the initial condition for rjt is instead of X[-n,n]<^- Indeed, using 

(2.2) with D = Z'^ gives 

Urn p''l-"M'i (At ^0\/t>O)= Urn Urn P^'l-"^"!'* (^^ 0) 

n— >oo n— >oot— >oo 

= Urn Mm P'^z'* (At n [-n, n]'^ 7^ 0) 

n— ♦00 >oo 

= Urn V({(E, F):En [-n, nf 7^ 0}) 

n — >oo 

= I7({(i?,F):£;/0}). 

This last probability is 1 when rjt survives, so in this case given any e > 
we can choose a positive integer m such that 

(4.1) p'^i-m.mjd (^4 / Vt > 0) > 1 - e. 

Recall that in Proposition 2.1 we showed that the limit distributions of 
the processes started at Xz'* ^^'^ ^z<^ same. It is then reasonable 

to expect that the asymptotic behavior as t ^ cxo of the process started at 
X[-n,n]'i is similar to that of the process started at J^[__„.„]ti, at least for large 
enough n. This idea will allow us to derive the lemma from (4.1). 

Proof of Lemma 4.1. Let e > and choose m to be the positive 
integer obtained in (4.1). To extend this inequality to the process started at 
X[-n,n]'*' we will consider two copies of the process r]t and 7]t coupled using 
the graphical representation, with rjj started at i^[_m,m]<* Vt X[-n,n]<* 
for some large n> m. For simplicity, we will write Q{k) = [—k,k]'^. 

We want to obtain a space-time cone growing linearly in time such that 
Ut>o{^} ^ A^'^'™' is contained in that cone with high probability. To achieve 
this, we compare with a branching random walk Zt with branching 

rate 13/ {2d) and no deaths [i.e., each particle in Zt gives birth to a new 
particle at each neighbor at rate (3/ {2d), and multiple particles per site 
are allowed]. Let {pt{x,y)}x^y£Z'' be the transition probabilities of a simple 
random walk in Z"^ that moves to each neighbor at rate 13/ {2d) and let Ct 
be the set-valued process given by 

Ct = {x(^r^:Zt{x)>Q}. 

For DCZ'^, Zf and will denote the processes started with all sites in 
D occupied by one particle and no particles outside L). It is not hard to see 
that for any t>0 and any rc G Z*^, 

E(zf^(x)) = e'3*Pi(0,x) 

[see, e.g., the proof of Proposition 1.1.21 in Liggett (1999)]. Therefore, for 
any D C Z'^, 

E(|cf > nD^\)<J2 n4'\x)) = e^' pti^^^)- 
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From this we get that if A; > m and c > then 

(4.2) Ei\cf^^"'^nQik + cty\)<i2m + l)^e''' ^ ptiO,x). 

\\x\\cx>>k—in+ct 

Now if Xt is the one dimensional random walk starting at and moving to 
each neighbor at rate (3/ {2d), Chebyshev's inequality gives 

¥{\Xt\ >k-m + ct) = 2F{Xt -k + m-ct>0) 

The last equality can be obtained by seeing Xt as the difference between 
two independent Poisson random variables, each with mean {(3t)/{2d), and 
using the fact that the moment generating function of a Poisson random 
variable Y with mean A is E(e*^) = e'*'^'^''"^). Applying this bound to each 
coordinate of the d-dimensional walk, we get that 

Pt{Q,x)<df{\Xt\>k-m + ct) 

\\x\\ao>k—m+ct 

and then using (4.2), we deduce that c can be taken large enough so that 
jg^l^QM ^ ^ ^^)C|^ < 2d(2m + l)^e-('^-'")e-*. 

Observe that, by the definition of Zt, the process is dominated by 

^Q(m) ^ ^^^^ bound implies that 

/ poo \ poo 

e( / lAt'^^"-^ nQ{k + cty\dt] < E{\cf^'^'^ nQ{k + cty\)dt 

(4.3) ^ 

<2d(2m + l)V('=-™). 

We can use this inequality to estimate the probability that At C Q(/c + 1 + 
ct) for all t > 0. Observe that if x € At n Q{k + 1 + ctY, the particle at x 
survives at least until time t + 2/c with probability e~'^"^^^^^^'^ , and thus 
X AgD Q{k + csY for all s G [t + 1/c, t + 2/c] with at least that probability. 
We deduce that 



E'^of™) {^j^\^t r\Q{k + cty\ dt 



> p'^Q(-) [At nQ{k + l + ctf / for some t > 0)e-2"(^+'^)/^-. 

c 

Therefore, if we let 

Gi = {Aj CQ{k + l + ct) yt > 0}, 
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(where denotes the set of I's in the process rjj started at i^Q^m))! we can 
use this bound together with (4.3) to get 

P(G5) < 2cd{2m + l)d^2ail+S)/c^-ik-m)_ 

Choosing now k large enough yields 

(4.4) P(Gi)>l-e. 

Now take n > fc, T > 0, let (i - r)+ = (t - T) V 0, and cah G2 the event 
that on the space-time region Uf>o{^} ^ + "^(^ ~ ^)^) the environment 
for rj^ dominates the environment for r]l [with respect to the order (1.1)]: 

G2 = {B^ C b} on Q{n + c(t - r)+) Vt > 0}. 

We want this space-time region to contain the region defining Gi, so we let 
T={n-k-l)/c. 

Observe that, since we are coupling the processes using the canonical 
coupling given by the graphical representation, once the environment is equal 
for both process at a given site, it stays equal at that site from that time on. 
In particular, dominates Bj on Q{n) for all t > 0. For any other site, any 
symbol •_! or *_i leaves the environment equal for both process. Therefore, 

P(G^)< IP(no or *_i at X by time r+ (||x||oo -n)/c) 

j>n 

^ ^a{i+5){k+i)/cj2(2j + i)'^e~"(^+'^)^'/''. 

j>n 

By taking n large enough, we obtain 

(4.5) P(G2)>l-e. 
Finally, let 

Gs = {Al^0yt>O}. 

By (4.1), P(G3) >l-£. Observe that on the event Gi n G2 n G3, contains 
I's at all times with probability 1. Therefore, 

pX[-n.n]d (^At / vt > 0) > P(Gi n G2 n G3) 

>1-P(G5)-P(G1)-P(G5;) 
> 1 - 3e, 

and choosing e small enough we get the result. □ 
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In the following lemma, we combine and extend for our process the results 
in Liggett (1999) that lead to the block conditions. Consider the process iVt, 
for L > 0, where no births are allowed outside of {—L,L)'^. Define N^[L,T) 
to be the maximal number of space-time points in 

S+{L,T) = {{x,s) e {{L} X [0,L)''-^) X [0,T]:x G lA^} 

such that each pair of these points having the same spatial coordinate have 
their time coordinates at distance at least 1. 

Lemma 4.2. Suppose that the process survives. Then for any a > 0, there 
is a positive integer n satisfying the following: For any given pair of positive 
integers N and M , there are choices of a positive integer L and a positive 
real number T such that 

(4.6a) P^[-"."l'' (Il^t n [0, Lf\ >N)>1- a^'" 

and 

(4.6b) pX[-„,„id ( jv+ (L, T) > M) > 1 - a'^'"!'^. 

Proof. By Lemma 4.1, we can choose a large enough integer n such 
that 

(4.7) P^[-"."l'*(^i/0 Vt>0)>l-o-^. 

Having this, the proof of the lemma is a simple adaptation of the corre- 
sponding proofs for the contact process. To avoid repetition of published 
results, we will explain the main ideas involved and why the original proofs 
still work with the random environment, but refer the reader to Section 1.2 
of Liggett (1999) for the details. 

We claim the following: For any finite A'^'L'^ and any > 1, 

(4.8) lim lim P'^-^di^il >iV)=P>^-4(-^^-^0 Vt>0). 

t — >oo L — ^oo 

To see that this is true, we observe that 

lim P^^-^d^Atl >iV) =P^-4(|^j| >N) 

and then argue that, conditioned on survival, \At \ — > oo as t — > oo with 
probability 1. This follows from the easy fact that there is an en > such 
that if 1^1 < A'^ then the process started with I's at ^ becomes extinct with 
probability at least e^r, so 

P'^-^(0 < \At\ < N)eN <F^^{t<T <oo) — > 0. 

t— »oo 

The next step is to use positive correlations to localize estimates on the 
cardinality of ^At to a specific orthant of Z'^: For every > 1 and L>n, 

(4.9) P^l-^.^l^di^tn [0,L)''| < Af) < [P^[-"."l''(|Ly4d <2'^A^)]^"''. 
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This relation follows easily from the positive correlations result in Lemma 
2.3, and its proof is identical the proof of Proposition 1.2.6 in Liggett (1999). 

Observe that (4.7), (4.8) and (4.9) together suffice to obtain (4.6a). The 
preceding arguments can be modified to obtain similar estimates for N^{L, T), 
which in turn give (4.6b). The only detail remaining is getting the same L 
and T to work for both inequalities. This is done by obtaining sequences 
Lj y oo and Tj y oo such that (4.6a) holds with L = Lj and T = Tj for 
every j > 1, and then adapting the arguments above to show that (4.6b) 
must hold for some pair (Lj ,Tj). We refer the reader to the proof of The- 
orem L2.12 in Liggett's book for the details on how this is achieved, and 
remark that the argument depends only on properties such as positive cor- 
relations and the Feller property which are available both for r]t and the 
contact process. □ 

We state now the block conditions that are equivalent to the survival of 
the process. 

Theorem 4.3. The process survives if and only if for any given e > 
there are positive integers n and L and a positive real number T such that 
the following conditions (BC) are satisfied: 

(BCl) P^[-"."l'' (i+2n,^T+i + [-n, n]'^ for some x £ [0, L)'^) > 1 - e 
and 

pX[_„,„j<i (^_^2n^t+i + [-n, nf for some 0<t<T 

(BC2) 

and some x £ {L + n} x [0, L) ) > 1 — e. 

Observe that these conditions correspond exactly to the conditions in 
Theorem L2.12 of Liggett (1999). This will allow us to borrow the argu- 
ments from Liggett's book to prove that (BC) implies survival for rjt. The 
reason why we need the conditions (BC) starting r]t from X[-n,n]<^ is because 
the proof of their sufficiency for survival (as well as their use in the proof 
of Theorem 2) demands obtaining repeatedly cubes fully occupied by I's 
and, at each step, restarting the process at the lowest possible configuration 
having those cubes fully occupied. 

Proof of Theorem 4.3. The proof uses the exact same arguments 
as those in the proofs of Theorems 1.2.12 and 1.2.23 in Liggett (1999). As 
before, we will only make some remarks and refer the reader to Liggett's 
book for the details. 

The necessity of (BC) follows from Lemma 4.2, by choosing the quantities 
N and M to be large enough to produce the desired boxes filled with I's. 
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For the sufficiency of (BC), attractiveness and (S2) imply that it is enough 
to show that for some n > the process started at X[-n,n]'' contains I's at 
all times (by using, as above, the fact that for any given n > the process 
started at X{o} [— n,n]'^ fully occupied by time 1 with some positive 
probability). The proof of this fact relies on starting with a large enough 
cube fully occupied by I's and then moving its center in an appropriate 
way. This is used to compare the process with supercritical oriented site 
percolation, and conclude that such boxes exist for all times with positive 
probability. □ 

The following consequence of Theorem 4.3 is obtained in the same way as 
for the contact process, see Theorem 1.2.25 in Liggett (1999) for the details. 



Corollary 4.4. If [5 = (3c{a,5) or 6 = 6c{a, f3), then the process dies 
out. 

5. Complete convergence. We are ready now to use the block construc- 
tion of Section 4 to prove Theorem 2. The key step in the proof will be 
to obtain the result in the special case where the initial distribution // is a 
probability measure of the form i'a, in which case we can use duality. 

Proposition 5.1. For every A c Z"', 

'q'^A ^ pi^A < J, _^ pi^A = 

To prove the proposition, we need a preliminary lemma. Both the proof 
of the proposition and this lemma are inspired by the proof Theorem 2 in 
Durrett and M0ller (1991). 

We will denote by W'^A''^c ^j^e probability measure associated to starting 
the process at i/a and its dual at i^c, using the same realization of the 
graphical representation, as explained in Section 2.2. 

Lemma 5.2. For every finite C C Z*^ and every e > 0, if r is a positive 
real number and s is large enough, then 



s 



< e. 



Observe that for the (ordinary) contact process, the forward process and 
the dual are independent when they run on nonoverlapping time intervals, 
so this fact is trivial and holds with s/2 replaced by s. 
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Proof of Lemma 5.2. Given r and e, there is a g = (?(|C|) such that 
every dual active path in (?}^'^''^)o<ti<r stays inside C + [—q,q\'^ with proba- 
bihty at least 1 — e. To see this, observe that the number of particles in all 
such dual active paths is dominated by Xr, where {Xr)r>o is a branching 
process starting with |C| particles and with birth rate P and death rate 
(we are ignoring deaths and coalescence of paths). By Markov's inequality, 
¥(Xr > q) <: E{Xr)/q < e for large enough q. Since any dual active path in 
fj^^'^ starts inside C, Xr < q implies that all dual active paths are contained 
inside C + [—q,q]'^ up to time r. 

Now denote by jy^^'"''/^^ and fj^^'"^^'^'^'^ modifications of the process and its 
dual, constructed on the same graphical representation as the original ones, 
where the environment is reset at time s/2 to its equilibrium /Xp, indepen- 
dently of its state before s/2 (that is, at time s/2 we replace every —1 by a 
and then flip every site to —1 with probability p, regardless of it being at 
state or 1). Then for given r and q, if s is large enough, we have that 

P'^-4''^c'(5^ = 5(mp,^/2) on C7+[-g,g]'^ VnG [s,s + r]) 

(5.1) 

> (1 _ ^-o.{l+5)s/2^\C+[~q,qY\ > 1 _ 

Indeed, for any given x G C + [—q^q]*^ the probability that and B^'"^^'^^ 
ar equal at x for every u G [s, s + r] is bounded below by the probability that 
an exponential random variable with parameter a{l + 6) is smaller than s/2 
[because any symbol or *_i above {x,s/2) leaves the environment at 
that site equal for both processes from that time on]. 

The property discussed at the first paragraph of the proof together with 
(5.1) imply that 



2 7 — /-T — II 2s. 

The statement of the lemma follows now from the independence of disjoint 
parts of the graphical representation and the stationarity of Bt , since 

^yA,yc > 1^ ^(Mp,s/2),r+s _^ 0^ ^(A.p,./2),r+. ^ ^ _^ ^ 

^^>0r^(i;^0,s5nz?/0). □ 



Proof of Proposition 5.1. The result is straightforward in the sub- 
critical case. If the process survives, and since weak convergence in this 
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setting corresponds to the convergence of the finite-dimensional distribu- 
tions, it is enough to prove that the following three properties hold for any 
two finite subsets C, D oil/: 

(cl) P''^(Aj nC/0) — > P'^-4(r = oo)F({(^,F):£;nC/0}), 

t— ►oo 

P'^'H^f n L» / 0) = P^^(r<oo)i/({(^,F):FnD/0}) 

(c2) 

+ P^^ (r = oo)v{{{E, F):FnD^0}) 

and 
(c3) 

— > F"^ (r = oo)F({(^, F):EnC^0,FnD^0}). 

►oo 

Indeed, all the finite-dimensional distributions of the process are determined 
by these probabilities via the inclusion-exclusion formula. Observe that the 
right-hand side of (c2) is equal to iJ,p{{r] -.rj^x) = —1 for some x G D}). 

The convergence in (cl) follows from the same arguments used in Liggett 
(1999) for the contact process. Using duality (Proposition 2.2), the proof of 
Theorem 1.1.12 in that book applies in the same way to obtain the fact that 
(cl) holds if and only if for every rr G Z*^ and every A C Z*^, 

(5.2a) P'"-4(r = oo)=P'''^(2;G^i i.o.) 

and 

(5.2b) lim liminf P''[-"."l'* (^ n [-n, n]'^ / 0) = 1. 

n— >oo t— *oo 

The analogous conditions are checked for the contact process in the proof 
of Theorem 1.2.27 in Liggett (1999). The equality in (5.2a) follows from the 
same proof after some minor modifications, so we will skip the argument. For 
(5.2b), Theorem 4.3 allows us to use Liggett's arguments to get the desired 
limit when i^[_„^„]<i is replaced by X[-n,n]d-^ so given any e > we can choose 
a large enough integer m such that 

(5.3) liminf P'^^I-™.-]'' {At n [-m, ^ 0) > 1 - e. 

t— >oo 

Given this m, we can choose a large enough n so that the process started at 
^l-n,n]'i contains at time a fully occupied cube of side 2m + 1 (contained in 
[— n,n]"') with probability at least 1 — e (in fact, any translate of [— m,m]°' 
contained in [—n,n\'^ is fully occupied by I's with some probability p > 0, so 
we only need to choose n so that [—n,n\'^ contains enough disjoint translates 
of [— m,m]'^). On this event, we can restart the process by putting every site 
outside that cube at state —1 and use attractiveness, translation invariance 
and (5.3) to get 

liminf P''l-"'"l''(At n [-n, n]'^ ^ 0) > {1 - ef, 

t— >oo 



24 D. REMENIK 

whence (5.2b) follows. There is only one detail to consider: in his book, 
Liggett only proves the condition analogous to (5.2b) in case d > 2, be- 
cause it is simpler and the case d=l was already done in Liggett (1985). 
The difficulty in the one-dimensional case arises from the fact that certain 
block events are not independent. This can be overcome by comparing with 
/c-dependent oriented site percolation instead of ordinary oriented site per- 
colation [see Theorem B26 in Liggett (1999)]. We refer the reader to Section 
5 of Durrett and Schonmann (1987), where the authors use a similar block 
construction to derive a complete convergence theorem for a general class of 
one-dimensional growth models. 

The convergence in (c2) is trivial due to the stationarity of the environ- 
ment process. To prove (c3), we start by observing that 

F-"^ {Ar+s n C / 0, Br+s r\D^0) 

(5.4) 

which follows from constructing (r/^'^)o<u<r+s and (f/^*^ )o<n<r+s on the 
same copy of the graphical representation. On the other hand, 

- P^^''^^ {As^0, AI+' / 0, Bl+' nD^0)\ 
(5.5) = P'^--*''^^ {As / 0, / 0, A, n = 0, B^+' n D / 0) 

< ^uA,uc {As^0, / 0, A, n AI+' = 0) 

= F^^^^^iAs / 0,i;+" / 0) - ^''^'''^{As n ^ 0). 

Observe that 

P'''^(s/2<r<oo) — ^0. 

s— >oo 

Thus, for any given D CZ,'^ and e > 0, and for large enough s, we can write 

IF^^'^c (^As / 0, ^ 0, nD^0) 

- F"^'"^ (r > s /2, A;+' / 0, ^5+^ nD^0)\ 

(5-6) 

= P^^'^c (s/2 < r < s, ^;+'' / 0, n D / 0) 
<P''-4(s/2<r<oo) < ^. 
Putting the previous observations together we get, for large enough s 

|P^^ {Ar+s n C / 0, Br+s nD^0) 
- F'^^'^'^iAs / 0,a;+' / 0,^^'+' n D / 0)1 
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= |p'^^''^c:(^^ n + 0, n L> 0) 

-P'^-^'^^(A, /0,i;+V 0,^5+^ nz)/ 0)1 by (5.4) 

< £ + \^^A.vc(j > s/2,i;+^ / 0) - P''-^'''^(A, n / 0)1 by (5.6), 

o 

where we used D = U^ and the fact that W-'^^''^ [b^^^ ^ 0) = l in the apph- 
cation of (5.6). Using again this fact to apply Lemma 5.2 with D = Z,"^, and 
then using duahty we get 



r > s/2, A;+' / 0) - F^^'^^iAs n / 0) 



< - + |P'^^(r > s/2)P'^c'(ir _^ J2;) - P'^^'^c(^^ p ^r+s ^ 
3 

= £ + > s/2)P''z'^ (A^ n c 0) - P^^(A,.+, n c / 0)1 

o 

for large enough s. By (cl), the last difference converges to as r, s — > 00, 
so we finally get 

IP'^^ {Ar+s n C / 0, Br+s nD^0) 
- F''^'''^ {As + 0, 0, ^^'+" n L> / 0) 1 < £ 

for large enough r, s. 

This calculation implies that in order to prove (c3) it is enough to show 
that 

^vA.vc ^ 0, ^ 0, ^5+^ n D / 0) 

— . P''^(r = oo)F((E,F):^nC/0,FnZ:'^0). 

Repeating the previous application of (5.6) and Lemma 5.2 we get that, for 
large enough s, 

\^^A,uc^J^^ ^ 0^ Xr+s ^ 0^ j^r+s n 2) / 0) 



F""^ (r > s /2)P^c (i; / 0, ^5 n i:* y^ 0) 
< £ + ip-^A.^c > ^ 0, n z? / 0) 

^'^■^(t > s/2)P''c (i;: / 0, ^5 n / 0) 



Therefore, we can finally reduce to proving that 

p^^ (r > s /2)P'^c' (i; / 0, ^5 n D / 0) 

^(r = oo)v{{E, F):EnC ^0,FnD^0). 
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This follows easily from duality, since (2.1) yields 
P^^(r > s/2)lP''^(i; / 0, n L> / 0) 

= P'^^ (r > s/2)P'"z'' {Ar n C / 0, s,, n / 0), 
and this last term converges to the desired limit as r, s — > oo. □ 

We extend now Proposition 5.1 to the general case. 

Proof of Theorem 2. It is enough to show that 

(5.7) hmE''(/(r?t))=P^(r<oo) J /dz^ + P^(r = oo) J f dV 

for every / in the space of continuous increasing functions depending on 
finitely many coordinates of X, which we will denote by J-. To see this, 
observe that given any two finite subsets C, D of Z*^, the functions 

fiiE,F) = Iedc^z, 

f2{E,F) = lFnD=0 and faiE, F) = lEnCyt0,FnD=0 

are all in and (as in the proof of Proposition 5.1) all the finite-dimensional 
distributions of the process can be obtained from E^(/i(?7j)), E'^(/2(r/t)), and 
&{f^{r]t)) by the inclusion-exclusion formula. 

Let / be a function in and observe that, in particular, / is bounded. 
One inequality in (5.7) is easy: by the Markov property and attractiveness, 
given < s < t we have that 

E'^(/(7yt)) = &if{7]t),T <s)+ W{f{r,t),T > s) 

= E^^iE^^ {f{rit^s)),r < s) + E^(E''= {f{vt^s)),r > s) 

< E(/(??o„J)P^(t < s) +E>^^'*(/(??t„,))P^(T > s), 

where 77° denotes the process started at the configuration r] = 0. Since rj^ =^ 
Hp = E and ?/f =^ V, we get 

limsupE''(/(^t)) <IP^(-^ < -5) / fdii + F^'{T>s) [ fdU, 
and now taking s — > 00 we deduce that 

(5.8) limsupE^(/(77j))<P''(r<oo) f fdu + F^{T = 00) f fdv. 

t — >oo J J 

To obtain the other inequality in (5.7), we will begin by considering the 
case // = X[-n,n]'* and showing that, given any e > and any x € Z'^, 

(5.9) liminfE^-+[-"."l''(/(r/t),T = oo) > [ f dV - e 
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for large enough n. By the translation invariance of r}t and 17, it is enough 
to consider the case x = 0. To show (5.9), we will use the construction in- 
troduced in the proof of Lemma 4.1. Using the notation of that proof, recall 
that we showed that, given any 7 > 0, there are positive integers n> k> m 
such that 

P(GinG2nG3)>l-37. 

This means that the processes rj^ (started at z^[_m,m]'*) ^^'^ Vt (started at 
X[-n,n]'') can be coupled in such a way that, with probability at least 1 — 87, 
for ah t > we have that Aj / 0, / 0, Aj C Q{k + 1 + ct) and Bf C B} 
inside Q{k + 1 + ct). 

Let G = Gi n G2 n G3 and 7 > and choose n > k > m so that P(G) > 
1 — 87. We will denote by and the extinction times of the processes rjl 
and r]^, respectively. Define 

Vt = {AlB!) 

and observe that, on the event G, fjt defines an ^-valued process and, more- 
over, 7]^ > fjt for all t>0. Therefore, since / is increasing and {r^ = 00} ^ G, 

(5.10) E(/(%'),r2 = oo)>E(/(7?2),G)>E(/(f?0,G) 
for all t>0. Now observe that, trivially, 

nf{flt),G) = nf{m),r^ = 00) - E(/(fyi), r2 = 00, G^), 

and 

nf{flt),r^ = oo,G^) < ||/||ocIP(G^) < 37II/II00, 

so 

(5.11) E(/(r?t), G) > E(/(r?t), = c^) - 87II/II00 
for all i > 0. On the other hand, 

(5.12) |E(/(r?t), r2 = 00) - E(/(r?,i), = oo)| 0. 

To see this, observe that since / depends on finitely many coordinates, then 
given any g > 0, f ijjs) = f {fll) for all s > t with probability at least 1 — g if 
t is large enough. Indeed, if C Z*^ is the finite set of coordinates of X on 
which / depends, then repeating the calculations that led to (4.5) we get 
that 

¥{B\{x) 7^ Bg{x) for some x & K and some s>t) 

< y P(no or *_i at x by time t) = |/C|e-°(^+'^)* — > 0. 

— ' i— »oo 
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Therefore, given any q> 0, 

mf{Vt),r^ = oo) - E(/(r?i),r2 = oo)| < E{\f{fit) - f{vl)\) < 2q\\f\\^ 

for large enough t, and we get (5.12). Finally, we have that 

E(/(7?i),r2 = oo) =E(/(r?i),ri = oo) - (E(/(7?i), = oo) - E(/(r?i), G)) 
-(E(/(r?i),G)-E(/(7?i),r2 = oo), 

and since G C {t^ = oo} n {r^ = oo}, 

mfivhy = oo) - E{f{vj),G)\ < ii/iiooip(G^) < 3711/1100 

for i = l,2. Thus, Proposition 5.1 implies that 

liminfE(/(7?i),r2 = oo)>P(r^ = oo) / /(fF - 67II/II00, 
and since F{t^ = 00) > P(G) > 1 — 87, we obtain 

(5.13) liminfE(/(r/i),r2 = oo)> f f dU - 9j\\f\\oo. 
Putting (5.10), (5.11), (5.12) and (5.13) together, we deduce that 

liminfE(/(7/2),r2 = oo)> / /(fI7 - 127||/||oo, 

and choosing 7 appropriately we obtain (5.9). 

Getting back to the proof of the remaining inequality in (5.7), let e > 
and choose n S N so that (5.9) holds. Define 

= inf{A; E N : r/^ ^ x + [— n, n]^ for some x G Z*^} 

and let p = (^Ai D x + [— n, n]'^ for some x € Z"^) > 0. Observe that for 
any k>0, if then A^+i contains some translate of [— n,n]'^ with 

probability at least p (by attractiveness and translation invariance) and, 
therefore, since the Poisson processes used in the graphical representation 
for disjoint time intervals are independent, we deduce that 

(5.14) {r = 00} C {iV < 00}. 

When N < 00, we will denote by X the center of the corresponding fully 
occupied box. If there is more than one point x such that x + [—n,n]'^ is 
fully occupied by I's at time A^, we pick X to be the one minimizing (j){x), 
where cf) is any fixed bijection between Z"^ and N (this ensures that the 
events {X = x} are disjoint for different x). Then given m G N, the Markov 
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property and attractiveness imply that 

m 

^''{firh),r = oo) > ^ &ifivt),r = oo, N = k) 

k=0 

m 

= ^ E^(E^^ {f{m~k).r = oc),N = k) 

k=0 

m 

>Y.Y. lE^(IE''-+[-"'"l'(/(^t-fc),r = oo),Af = A;,X = x) 

for t > m. Since / is bounded, (5.9) implies that 

liminfE'^(/(r/t),r = oo) >( f f dU - e) £ ^ P^(iV = A:, X = x) 

= (^j fdi7-e^F''{N<m). 
Taking now m — > cxd, we get by (5.14) that 

liminfE'^(/(r?t),r = oo) > (^J fdu-e^¥^{N< 



oo 



> 



if e < / / dU, and taking e ^ we deduce that 

liminfE''(/(r?t),r = oo) >P''(r = oo) / fdv. 

i-»oo J 

On the other hand, by arguments similar to those that led to (5.8) (using 
attractiveness to compare with the process started at X0)i we get 

liminfE^(/(??i),r <oo) >P^'(^<oo) / f di^- 

We finally deduce that 

liminfE^(/(r/t)) >P''('^< oo) J /(fi7 + P^(r = oo) J fdU, 

and the proof is ready. □ 

Acknowledgments. I am grateful to my advisor, Rick Durrett, for his 
guidance during the work in this paper. I thank him especially for his careful 
reading of several versions of the manuscript and the thorough comments 
and suggestions he provided, which led to a substantial improvement of the 
paper. I also wish to thank Ted Cox for reviewing part of an earlier version 
of the manuscript. Finally, I thank an anonymous referee for many detailed 
comments and suggestions which greatly improved the exposition of the 
paper. 



30 



D. REMENIK 



REFERENCES 

Andjel, E. D. (1992). Survival ol multidimensional contact process in random environ- 
ments. Bol. Soc. Brasil. Mat. (N.S.) 23 109-119. MR1203176 

Athreya, K. B. and Ney, P. E. (1972). Branching Processes. Die Grundlehren der 
mathematischen Wissenschaften 196. Springer, New York. MR0373040 

Bezuidenhout, C. and Grimmett, G. (1990). The critical contact process dies out. Ann. 
Probab. 18 1462-1482. MR1071804 

Bramson, M., Durrett, R. and Schonmann, R. H. (1991). The contact process in a 
random environment. Ann. Probab. 19 960-983. MR1112403 

Broman, E. I. (2007). Stochastic domination for a hidden Markov chain with applications 
to the contact process in a randomly evolving environment. Ann. Probab. 35 2263-2293. 
MR2353388 

Durrett, R. and M0ller, A. M. (1991). Complete convergence theorem for a compe- 
tition model. Probab. Theory Related Fields 88 121-136. MR1094080 

Durrett, R. and Schinazi, R. (1993). Asymptotic critical value for a competition model. 
Ann. Appl. Probab. 3 1047-1066. MR1241034 

Durrett, R. and Schonmann, R. H. (1987). Stochastic growth models. In Percolation 
Theory and Ergodic Theory of Infinite Particle Systems (Minneapolis, Minn., 1984- 
1985). The IMA Volumes in Mathematics and Its Applications 8 85-119. Springer, New 
York. MR894544 

Durrett, R. and Swindle, G. (1991). Are there bushes in a forest? Stochastic Process. 

Appl. 37 19-31. MR1091691 
Harris, T. E. (1972). Nearest-neighbor Markov interaction processes on multidimensional 

lattices. Adv. Math. 9 66-89. MR0307392 
Harris, T. E. (1974). Contact interactions on a lattice. Ann. Probab. 2 969-988. 

MR0356292 

Klein, A. (1994). Extinction of contact and percolation processes in a random environ- 
ment. Ann. Probab. 22 1227-1251. MR1303643 

Liggett, T. M. (1985). Interacting Particle Systems. Grundlehren der Mathematischen 
Wissenschaften [Fundamental Principles of Mathematical Sciences] 276. Springer, New 
York. MR776231 

Liggett, T. M. (1992). The survival of one-dimensional contact processes in random 
environments. Ann. Probab. 20 696-723. MRl 159569 

Liggett, T. M. (1999). Stochastic Interacting Systems: Contact, Voter and Exclusion 
Processes. Grundlehren der Mathematischen Wissenschaften [Fundamental Principles 
of Mathematical Sciences] 324 Springer, Berlin. MR1717346 

Luo, X. (1992). The Richardson model in a random environment. Stochastic Process. 
Appl. 42 283-289. MR1176502 

Meester, R. and Roy, R. (1996). Continuum Percolation. Cambridge Tracts in Mathe- 
matics 119. Cambridge Univ. Press, Cambridge. MR1409145 

Center for Applied Mathematics 

Cornell University 

657 Rhodes Hall 

Ithaca, New York 14853 

USA 

E-mail: dir4(acornell.edu 



